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Abstract 

^vq I We discuss the choice of input parameters for the renormahzation of the chargino 

and neutrahno sector in the minimal supersymmetric standard model (MSSM) in 

the on-shell scheme. We show that one should chose the masses of a bino-like, 

a wino-like and a higgsino~like state as inputs in order to avoid large corrections 

^^ I to the masses of the other eigenstates in this sector. We also show that schemes 

where the higgsino-like input state is a neutrahno are more stable than those where 

' r~|', the mass of the higgsino~like chargino is used as input. The most stable scheme 

Qh| uses the masses of the wino-like chargino as well as the masses of the bino- and 

p i' higgsino-like neutralinos as inputs. 

> 

oo 
(N 

o 



X 



\o 



*Authors' names are listed in alphabetical order 

tOn leave from Institute of Nuclear Physics, PAN, Krakow, ul. Radzikowskiego 152, Poland. 

1 



1 Introduction 

Among possible solutions of the hierarchy problem [H [21 [Sj H] , supersymmetric extensions 
of the Standard Model (SM) of particle physics [5], E] have the advantage that effects of 
the hypothetical new "superparticles" required in these constructions can be computed 
perturbatively. This distinguishes them e.g. from scenarios where electroweak symmetry 
breaking is accomplished by some new strong dynamicsHI At the same time, future 
measurements at the LHC [7J and a possible e^e~ linear collider [8] can easily have 
statistical uncertainties at or below the percent level, making correspondingly accurate 
theoretical calculations highly desirable. This will require the calculation of electroweak 
corrections at least at the one-loop level. Independent motivation for such calculations 
comes from the possibility that neutralinos form the dark matter in the universe, since the 
determination of the overall dark matter density of the universe has now been determined 
with an error of a few percent [H]. 

In all these cases the renormalization of the chargino and neutralino sector plays a 
central role. In most SUSY breaking scenarios the lighter charginos and neutralinos are 
among the lighter superparticles, i.e. among the first ones to be discovered by future e^e~ 
supercolliders. At the LHC the most important search channels start from the production 
of squarks and gluinos, which subsequently decay into charginos and/or neutralinos. 

The renormalization scheme of choice is on-shell renormalization, for (at least) two 
reasons. First, there is a fairly direct relation between the physical chargino and neu- 
tralino masses and the parameters that need to be renormalized; this relation becomes 
trivial in the limit where the spontaneous breaking of the electroweak gauge symmetry 
can be neglected. Moreover, masses will probably be among the first quantities that are 
measured accurately. 

The on-shell renormalization of the neutralino and chargino sectors in the minimal 
supersymmetric standard model (MSSM) with conserved i?— parity has been discussed 
in [TOl [TT| IT2] . In full generality, the chargino and neutralino mass matrices depend 
on six independent parameters that need to be renormalized. However, two of these 
(the masses of the W and Z bosons) already appear in the SM; they are renormalized 
through the same conditions as in the SM (albeit with additional contributions to the 
relevant two-point functions with superparticles in the loop). A third parameter, the 
ratio of vacuum expectation values (VEVs) tan/3, is usually renormalized in the MSSM 
Higgs sector. This leaves only three independent parameters whose renormalization has 
to be fixed through on-shell conditions for charginos and neutralinos. This means that 
only three of the six chargino and neutralino masses can be used as independent input 
parameters; their values are exact to all orders in the on-shell scheme. The other three 
masses can then be predicted, but the numerical values of these predictions are subject 
to loop corrections. 

Previously different calculations used different sets of input masses, corresponding to 
different variants of the on-shell scheme. For example, in ref . (TU] , both chargino masses 
and the lightest neutralino mass are used as inputs. In [13] the decay of the next to 
lightest neutralino to the lightest neutralino was studied at one loop level using the on- 
shell renormalization scheme. The input masses were conveniently assumed to be the 



•l-The "hidden sector" required to break supersymmetry may have non-perturbative couplings, but, 
being hidden, this generaUy does not invahdate perturbative calculations in the visible sector. 



masses of the two lightest neutrahnos and the heaviest chargino. 

Since the hghter neutrahno states are more accessible at colliders, using their masses 
as inputs seems favorable from the experimental point of view. However, reliable per- 
turbative predictions can be made only if the perturbative expansion is stable, i.e. if 
loop corrections are not very large. From this point of view using the masses of both 
charginos and the lightest neutralino as input can be problematic: the corrections to the 
masses of some of the other three neutrahnos can become large if the lightest neutralino 
is an almost pure higgsino or wino [HI [15]. The reason is that all three input masses 
then only depend very weakly on the U{1)y gaugino (bino) mass parameter. A large 
(finite part of the) counterterm may therefore be needed in order to correct small loop 
corrections to the input masses, so that these masses remain unchanged, as required in 
on-shell renormalization; the large counterterm will then lead to large changes of other 
neutralino masses. By a very similar argument, in this situation small experimental er- 
rors on the input masses can lead to very large errors on the bino mass parameter even 
when quantum corrections are neglected. An analogous instability can arise when using 
the masses of the two lightest neutrahnos and the heavier chargino as inputs. Note also 
that both these schemes require the mass of the heavier chargino as input, which is one 
of the heaviest states in this sector; this considerably weakens the argument in favor of 
using only the masses of the lighter neutralino states as inputs. 

If the mass of the higgsino-like chargino is used as input an additional perturbative 
instability occurs if the SU{2) gaugino (wino) and higgsino mass parameters are of similar 
magnitude. It has been noticed in fiU\ ITU [15] that using both chargino masses as input 
leads to undefined renormalized masses for exact equality of the absolute values of the 
wino and higgsino mass parameters; not surprisingly, corrections remain very large, albeit 
finite, if the difference between these two parameters is nonzero but small. We find that 
this instability also occurs in schemes where the mass of the higgsino-like chargino and 
two neutralino masses are used as inputs. 

In this paper we address these two stability issues by explicitly calculating chargino 
and neutralino masses at one-loop level using the on-shell renormalization scheme. We 
compare seven different choices for the three input masses, combining two, one or zero 
chargino masses with either the mass(es) of the lightest neutralino state(s) or with the 
masses of neutralino states of given character (bino-, wino- or higgsino-like). We demon- 
strate explicitly that the choice of any three input masses, which include a higgsino-like, 
a wino-like and a bino-like mass eigenstate, avoids the first instability, independent of 
whether chargino or neutralino masses are used. Requiring in addition perturbative sta- 
bility also for similar magnitudes of the wino and higgsino mass parameters singles out 
schemes where the masses of the (most) bino- and higgsino-like neutrahnos are used as 
inputs. Schemes using three neutralino masses as inputs show an additional, but much 
weaker, instability for approximately equal SU{2) and f/(l)y gaugino masses; this second 
instability is absent if the mass of the wino-like chargino, rather than that of the wino-like 
neutralino, is used as input together with the two neutralino masses listed above. 

In this paper we assume all relevant mass parameters to be real. Experimental con- 
straints on CP violation, mainly from EDM measurements, require their phases to be 
small, unless superpartners of first (and second) generation fermions are very heavy 

[lain]. 

This article is organized as follows. In Section [2] we discuss briefly seven different 



variants of the on-shell scheme, in the context of the MSSM with exact i?— parity. Sec- 
tion [3] discusses the practical implementation of our schemes. In SectionlHwe compare the 
above seven schemes, defined by different sets of input masses. Explicit expressions for 
the counterterms in schemes where two or three neutralino masses are used as inputs are 
given in Appendices A and B, respectively, while Appendix C shows how to distinguish 
the wino- and higgsino-like chargino state using Zx'^x'^ couplings. 

2 Formalism 

In this section we discuss the on-shell renormalization of the chargino and neutralino 
sector of the MSSM with exact i?— parity [lOl |T5]. To this end we first review tree-level 
results for chargino and neutralino masses and mixings. We then describe general features 
of the on-shell renormalization scheme as applied to this sector, before introducing three 
variants of this scheme, each of which has two or three subvariants. 

2.1 Tree— Level Results 

The tree level mass terms for the charginos, in the gauge eigenbasis, can be written as [5] 

- ^Lss = r^M^ij^ + h.c. (1) 

where ^^ ^^ 

,P+ = {W+, htf, ^-' = {W-, h^f (2) 



are column vectors whose components are Weyl spinorq^]. The mass matrix M*^ is given 
by 

V V2Mw COS 13 /i y ■ ^ ^ 

Here M2 is the SUSY breaking SU{2) gaugino (wino) mass, /i is the supersymmetric 
higgsino mass, M^ is the mass of the W boson, and tan /3 is the ratio of VEVs of the 
two neutral Higgs fields of the MSSM. M^ may be diagonalized using unitary matrices U 
and V to obtain the diagonal mass matrix, 

M' = U*M'V-' = I "l^^ ° 1 . (4) 

y m^+ J 

Without loss of generality, we order the eigenstates such that < m -+ < m-+ . Note that 

scenarios with m-+ ^ 100 GeV are excluded by chargino searches at LEP [18] . The left- 
and right-handed components of the corresponding Dirac mass eigenstates, the charginos 
xt with i = 1 or 2, are 

PLxi = v.j^l, p^xt = u:;^, (5) 

where P^ and Pr are chiral projectors, ^~ is the hermitean conjugate of the Weyl fermion 
ip7 , and summation over j is understood in Eqs.©. 



^Note that W~^ is the antiparticle of W , but /i^ is not related to h^ : the latter two fields reside in 
the two distinct Higgs superfields required in the MSSM. 



In the gauge eigenbasis, the tree level neutralino mass terms are given by [5] 

1 



-Cl^s = -r'M-r + h.c., (6) 



where 

T 



r=[BO, W^ hi hi) , (7) 

is again a column vector whose components are Weyl fermions, and the neutralino mass 
matrix M° is given by 



M° 



/Ml -Mzswcp Mzswsfi \ 

Ms Mzcwcp -Mzcwsp 

-Mzswcp MzCwC(s -fi 

\ MzSwS(3 -Mzcw-sp -/i / 



(8) 



Ml = N*M''N-^ 



m^o 
o' m^o 



Here sw,Si3,cw and c^ stand for sin ^vf, sin /9, cos 6*1^ and cos (3, respectively, where 9w is 
the weak mixing angle. Mz is the mass of the Z boson, and Mi is the SUSY breaking 
U{1)y gaugino (bino) mass. Finally, M2 and yU already appeared in the chargino mass 
matrix Q. 

M^ may be diagonalized by a unitary matrix A^ to obtain the diagonalized mass 
matrix M^, 

f m^o \ 

(9) 

V m^o / 

Again, without loss of generality, we order the eigenvalues such that 

< m^o < m^a < m^o < m^o. 

— Xi — X2 — X3 — X4 

Note that an arbitrarily light neutralino is still phenomenologically possible, if no relations 
between Mi, M2 and/or /i are imposed [19]. 

The left-handed components of the corresponding mass eigenstates, described by 
four-component Majorana neutralinos Xi with i G {1,2,3,4}, may be obtained as, 

PlX- = N^.i^l (10) 

where summation over j is again implied; the right-handed components of the neutralinos 
are determined by the Majorana condition xf = Xiy where the superscript C stands for 
charge conjugation. 

The parameters M2 and P can always be made real and positive through phase ro- 
tations of fields. If the parameters Mi and /i are also real, as we assume, the chargino 
mixing matrices U and V can be chosen to be purely real. However, the neutralino mixing 
matrix N can only chosen to be real if negative neutralino "masses" (better: eigenvalues 
of M°) are tolerated. Physical masses must not be negative, of course. In that case the 
entries of N can be either real or purely imaginary. Genuinely complex entries of N, as 
well as U and V, are only required if Mi and/or fi have non-trivial CP-odd phases. 



2.2 Generalities of On— Shell Renormalization 

At the quantum level, diagrams with internal and/or external chargino and/or neutralino 
lines will receive two-point function corrections on all these lines. These corrections are in 
general infinite. In order to absorb such infinities through renormalization, the mass ma- 
trices and the mass eigenstates need to be redefined. Here we discuss this renormalization 
at the one-loop level. 

The one-loop mass eigenstates can be related to the tree level mass eigenstates through 
wave function renormalization: 

X.'^^' = ('5., + \5Z.,Pl + \6Z:^Pr) ;e/-°™^i-<=^ (11) 

This relation holds for both the chargino sector, with Xi — xt i ^ ^ {l^^}, and the 
neutralino sector, with Xi = x°, i G {1,2,3,4}. Note that the wave function renormal- 
ization matrices 5Z'^ and 5Z" in the chargino and neutralino sectors are not diagonal. In 
the on-shell renormalization scheme loop-induced changes of the identities of the mass 
eigenstates can be described entirely through these wave function renormalizations. The 
entries of the mixing matrices N, U and V can therefore be chosen to be identical to 
their tree level values [15]. Explicit expressions for the wave function renormalization 
constants are not needed for the calculation of the chargino and neutralino masses, on 
which we focus in this paper. 

These masses receive explicit corrections from one-loop diagrams, which can be writ- 
ten as 



1 



1 



where / denotes a fermion species with mass rrif. 9^e denotes the real parts of the loop 
integrals involved, leaving imaginary parts of couplings unchanged (in the CP-conserving 
case these can only appear through the neutralino mixing matrix A^, as remarked at the 
end of the previous Subsection). Moreover, the S refer to various terms in the general 
two-point function of fermion / in momentum space: 

S//(P) = ^ [PlS^^(p) + P^S^^(p)] + Pl^'\p) + PrJ:'\p) , (13) 

where Pl and Pr are again the chiral projectors. Note that only diagonal two-point 
functions contribute to the corrections to physical masses at one-loop level in the on- 
shell scheme O 

The corrections of Eq. flT2|) are in general divergent. These divergencies are absorbed 
into counterterms to the mass matrices M^ and M"^ of Eqs.([3]) and ([8]), respectively, i.e. 

j^h^ve ^ prenormalized _^ ^j^_ ^^^^ 

As usual, we define the physical (on-shell) masses as poles of the real parts of the (one- 
loop corrected) propagators. The physical chargino masses are then given by, 

ml% = m-+ + {U*6M'V-% - 6m-+ ; (15) 



"OfF-diagonal two-point function corrections do appear in generic one-loop diagrams. Their infinities 
are absorbed in off-diagonal wave function renormalization constants, which are fixed by the definition 
that particles do not mix on-shell. See e.g. ref.[TOj for further details. 
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the corresponding expression for the neutrahnos is 

mg = m^o + {N*5M''N-^)ii - 6m^o . (16) 

The masses m-+,o appearing on the right-hand sides of Eqs. (IT5l) and (TT6l) are the (finite) 
tree-level masses. U, V and A^ are the mixing matrices in the chargino and neutralino 
sectors, respectively, which, as already noted, do not get modified by loop corrections. 
Sm-+.o are the explicit loop corrections of Eq. (IT^ as applied to the charginos and neu- 
tralinos. Finally, SM'^'^ are the counterterm matrices of Eq. flT^ for the chargino and 
neutralino sector, which we yet have to determine. 

Renormalizability requires that the counterterm matrices have exactly the same form 
as the tree-level mass matrices, i.e. all non-vanishing entries are replaced by the cor- 
responding counterterms while the vanishing entries receive no correction. To one-loop 
order counterterms to products like Mw sin /3 can be written as {6Mw) sin f3 + MwS sin /3, 
and so on. 

Altogether there are thus seven different counterterms: 6M]y, SMz, S9w, 5tan/3, SMi, 
SM2 and 6fi. The first three of these already appear in the SM. We renormalize them 
according to the on-shell prescription of electroweak renormalization, where Mw and Mz 
are physical (pole) masses, and cos 9w = Mw/Mz- This gives [20] : 

6M^ = mcEwwiM^); 
5Ml = ^t^zz{Ml); 

Mw fSMw SMz\ .^„. 

Here T,ww and T,zz are the transverse components of the diagonal W and Z two-point 
functions in momentum space, respectively; again only the real parts of the loop functions 
should be included as indicated by the 9^c symbols. Note that, while the definitions of 
these three counterterms are formally as in the SM, in the MSSM there are many new 
contributions to T^ww and T,zz involving loops of superparticles and additional Higgs 
bosons. 

The counterterms 01 7p imply that the sum of the squares of the VEVs of the two 
neutral Higgs components is already fixed. However, the ratio of these VEVs, tan/3, 
does not affect the W and Z masses, hence its counterterm has not yet been determined. 
Following refs. [211 ES] we fix it by the requirement that the two-point function connecting 
the CP-odd Higgs boson A to the (longitudinal part of the) Z— boson vanishes when A 
is on-shell. This giveqj 

aan/3 = ^^^^-l-^am(S^^(mi)). (18) 

Hence only the counterterms to Mi, M2 and /i remain to be fixed in the chargino and 
neutralino sector. This means that we can only require that three of the six physical 



*Note that the couphngs of A contain an extra factor of i relative to gauge couphngs. Therefore the 
imaginary part of the two-point function appears in Ea.(|18p: this contains the real (dispersive), infinite, 
part of the loop function. 



masses in this sector are not changed by loop corrections, i.e. only three of the six "tree- 
level" masses are physical (all-order) masses. The other three masses will receive finite, 
but non-zero correctionsJH In order to complete the definition of our renormalization 
scheme, we have to decide which three masses to use as input masses. There are many 
ways to do so; this is the topic of the following Subsection. 

2.3 Choosing Input Masses 

/ 6 \ 
For a fixed point in parameter space, there are =20 different ways to select three 



^ 3 ^ 

out of six masses. However, as we will see shortly, any scheme that always chooses masses 
with fixed subscripts i,j, k as inputs is bound to lead to perturbative instability in large 
regions of parameter space. Instead, one should choose the masses of chargino and/or 
neutralino states with specific properties as inputs; the indices of the input states will 
then vary over parameter space. 

We illustrate this by considering three different schemes, each of which has two or 
three variants. The most widely used inputs are the masses of both charginos and of 
the lightest neutralino fiU\ [15] . We call this scheme la. Using both chargino masses 
as inputs guarantees that the finite parts of the counterterms 6M2 and 5/i are usually 
small, since at least one chargino mass will have strong sensitivity to either of these 
counterterms. However, as already mentioned in the Introduction, this choice leads to 
instabilities if |A^ii| ^ 1 [HI [IS]- In this case the tree-level value of m^o depends only 
very weakly on Mi. Correspondingly the counterterm 6 Mi may require a very large 
finite part to cancel the finite part of the explicit corrections to m^o in Eq. fll6p . This in 
turn will lead to a very large correction to the mass of the (most) bino-like neutralino|i| 
In terms of parameter space, this instability will occur whenever |Mi| — M2 ^ Mz or 
I Ml I — |/x| ^ Mz, in which case Xi is wino- or higgsino-like. 

A suitable alternative is to instead use the mass of the bino-like neutralino as an input 
[T11[T5|. We call this scheme lb. In the limit of small mixing, this mass is approximately 
given by |Mi|, and it will retain strong sensitivity to Mi even in the presence of significant 
mixing. 

As already mentioned, in the on-shell scheme the input masses are interpreted as 
exact physical masses, i.e. their total one-loop corrections should vanish. Eqs. (ll5l) and 
([T6|) then lead to the following equations: 

{U*6M'V-^)n = Sm^+, 
iU*5M'V-%2 = Sm^p 
{N*6M''N-^)ii = 6m^o (no summation), (19) 

where i = 1 in scheme la, while i denotes the (most) bino-like neutralino in scheme 
lb. Note that both chargino equations in ([T9l) depend on both 6fi and 6M2, as well as 



TThe finiteness of these corrections is a condition for the renormahzabihty of the theory. In practice 
it affords non-trivial checks of our calculation. 

■I- Symbolically, if dm^^^f^/dMi — e, and explicit finite loop corrections to m^o are of typical (relative) 
one-loop order, i.e. O^a/n), then the finite part of the counterterm to Mi will be of relative order 
0{a/{eTT)), where a is an electroweak fine structure constant. The finite correction to the mass of the 
most bino-like neutralino will then also be enhanced by 1/e, and will thus become very large ii e^ a. 
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on 5Mw and (5tan/3 which have aheady been determined in Eqs. fll7l) and flTSl) . respec- 
tively. Similarly, the third (neutralino) equation in f lT9|) depends on all the counterterms. 
However, since the dependence on counterterms is linear, eqs. fllQI) can readily be solved 
analytically. 

Alternatively, one can use the masses of one chargino and two neutralinos inputs. The 
counterterms 5 Mi, SM2 and 5/i are then determined by the following equations: 

{U*5M^V-^)kk = (5m -+ (no summation). (20) 

In scheme 2a we take i = l,j = 2 and k = 1, i.e. the lightest chargino mass and 
the two lightest neutralino masses are used as inputs. These are the lightest particles in 
this sector, hence their masses are likely to be the first to be determined experimentally. 
However, this scheme is likely to lead to perturbative instabilities whenever the differences 
between |Mi|, M2 and |/i| are large. For example, if M2 — |/i| S> Mz, xt and X2 '^ill 
both be higgsino-like, and none of the input masses will be sensitive to M2, leading to 
a potentially large finite part of 6M2. If |/i| — M2 ^ Mz and |/i| — |Mi| ^ Mz, none 
of the input states is higgsino like; hence none of the input masses depends sensitively 
on /i, leading to a potentially large finite part of S/i. In scheme 2b i and j therefore 
denote the bino- and wino-like neutralino, respectively, while k denotes the higgsino- 
like chargino. This ensures that there is at least one input mass that is sensitive to 
each counterterm. This is true also in scheme 2c, where i and j denote the bino- and 
higgsino-like neutralino, respectively, while k denotes the wino-like chargino. 

With three neutralino masses as inputs, the counterterms 6M1, 6M2 and 6fi are all 
determined from Eq. flT6l) : 



{N*5M''N-^)ii = 5m^o, 
{N*6M'^N~^)jj = 6m^o, 
{N*5M''N-%k = (5m^o (no summation). (21) 

In scheme 3a we simply take the three lightest neutralinos as input states, i.e. i = 1, j = 
2 and k = 3 in Eqs. (l2T]) . This runs the risk that two of the input states are higgsino-like, 
in which case the third input mass cannot determine both 6M1 and 6M2 reliably. For 
example, if M2 — \Mi\ ^ Mz and M2 — |/u| ^ Mz, none of the input masses is sensitive 
to SM2, since the wino-like neutralino is the heaviest one. In scheme 3b i, j and k 
therefore denote the bino-like, wino-like and higgsino-like neutralino, respectively. 

For each set of input masses, the corresponding set of linear equations can be solved to 
obtain 6M1, 6M2 and 6fi. The solutions, in the context of scheme la, have been obtained 
using the mixing matrices in [TU] and without using the mixing matrices in [HI ITS] . The 
same solution can be used for scheme lb, with the substitution A^'iq, — )■ Nia, where i 
denotes the bino-like neutralino. Explicit solutions for schemes 2 and 3 can be found in 
Appendices A and B, respectively. 

The seven schemes are summarized in Table [H The second column lists the states 
whose masses are used as inputs. Here Xb, Xw and Xh stand for a bino-, wino- and 



Scheme 


Input states 


Regions of instability (counterterm) 


la 


Xi , X2 ' Xi 


Mi|-|/i|»M^(5Mi) 
Mi|-M2>Mz (5Mi) 


lb 


Xl > X2 ) X6 




2a 


Xi , Xi, X2 


M2- fi -:$> Mz {5M2) 
/i - M2 > Mz and |/x| - |Mi > Mz (5/i) 
Mi|-|/i|>Mz ((5Mi) 


2b 


Ah ' Xb? X«) 




2c 


Xui' Xb) X/i 




3a 


X?, xl xl 


M2 - /i > M^ and Mi| - /i > M^ (5Mi or 6M2) 
M2 - ;u > Mz and |Mi| < ^ {6M2) 
Mil - |/i| > Mz and M2 ^ /i (5Mi) 


3b 


~o ~o ~o 

Xb, X«,, X/j 





Table 1: Summary of the schemes discussed in this paper. The second column lists the 
states whose masses are used as inputs; here Xb, Xw and Xh stands for a bino-, wino- and 
higgsino-like state, respectively. The last column lists the regions of parameter space 
where schemes la, 2a and 3a may become perturbatively unstable, together with the 
counterterm that is poorly determined in that region. In the first instability region of 
scheme 3a, only one linear combination of 6M1 and 6M2 is well determined, depending 
on the relative ordering of |Mi| and M2. Note that the different regions of instability 
may overlap. 

higgsino-like state, respectively. The third column lists regions of parameter space, de- 
fined through strong inequalities, where at least one counterterm is poorly determined, 
leading to potentially large corrections to some mass(es); the poorly determined coun- 
terterm is given in parentheses. These regions only exist for schemes la, 2a and 3a; in 
fact, one can easily convince oneself that such regions of potential perturbative instability 
exist in all schemes that fix the indices of the input states a priori, independent of the 
characters of these states. 

This concludes the description of the renormalization schemes we are using. Before we 
turn to numerical results, we briefiy discuss some issues that might arise in the practical 
implementation of our formalism. 

3 Practical Considerations 

In this Section we discuss some issues regarding the practical implementation of our 
calculation. In particular, we define the bino-, wino- and higgsino-like states in terms of 
observable quantities. We then discuss how to extract the on-shell parameters Mi, M2 
and fi from the three input masses in the three schemes. Finally, we describe how to 
combine our formalism with spectrum calculators. 
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3.1 Definition of the States 

We saw at the end of the previous Section that schemes where the three states whose 
masses are taken as inputs include a bino-hke, a wino-hke and a higgsino-hke state are 
better behaved than schemes where the indices of the input states are fixed a priori. This 
raises the question how these input states are defined. 

We begin with the neutrahno sector. The bino-hke state is defined as the neutrahno x^ 
with the largest bino component, i.e. |A^ii| > l^jil Vj 7^ i. This can (at least in principle) 
be determined experimentally by finding the neutrahno with the strongest coupling to 
right-handed electrons, i.e. with largest absolute value of the Xi^r^r coupling]^ 

The (most) wino-like neutrahno is the state Xi* with largest wino component, i.e. 
1^42! > l^jil ^j 7^ i- Having already determined the bino-like neutrahno, the wino-like 
state is the remaining neutrahno with the strongest coupling to left-handed sleptons, i.e. 

the (non-bino) state with largest absolute values of the x^clGl and x^vl^l couplingsjj 

2 I I Ar \2\ 



Finally, we define the (most) higgsino-like state as the one with largest (|A^j3| +|A^, 



i4 



Unfortunately this sum does not directly correspond to a measurable coupling. However, 
using unitarity of the mixing matrix we can write jA^isp + |iVj4p = 1 — |-^iip — |-^i2p, 
i.e. the most higgsino-like state is the one with weakest couplings to first generation 
leptons. Moreover, once the bino- and wino-like states have been identified, in most 
cases it does not matter very much which of the two remaining states is defined to be 
"the" higgsino-like stateHIJ 

In the chargino sector, one can define the wino-like state as the state which has 
larger coupling to a (s)neutrino or to a left-handed (s)electron; the other state is then 
higgsino-like. Alternatively, one can take the strength of the xtx7 ^^ coupling as defining 
property: since the wino is an SU{2) triplet while the higgsino is a doublet, the former 
couples more strongly to the Z. We show in Appendix C that this also allows a unique 
definition of the wino- and higgsino-like states. 

3.2 Extraction of the On— Shell Parameters 

Experimentally one will (hopefully eventually) measure several neutrahno and chargino 
masses. In order to apply our formalism, one needs to calculate the input values of Mi, M2 
and n from the three input masses. This is necessary in order to predict the other three 
chargino and neutrahno masses, which were not used as inputs. This inversion from 
masses to parameters is also required when comparing results of different schemes (see 
below). 



^Occasionally there will be two neutralinos Xi'j X? with |A^ii| ~ |^ji| — l/v^- In this case it doesn't 
matter whether m^o or m^o is chosen as input, since both masses are then quite sensitive to Mi. 

Xi Xj 

^Occasionally the state x? with the largest \Nii\ also has the largest \Ni2\, i.e. the most bino-like 
state is also the most wino-like of all four states; this can happen in the presence of strong bino-wino 
mixing , i . e . Mi ~ M2 , if in addition 1 1 /i | — M2 \ ^ Mz ■ This ambiguity is resolved by first determining the 
most bino-like of all four neutralinos, and then determining the most wino-like of the remaining three 
neutralinos. A similar situation may be encountered in case of strong bino-higgsino mixing, especially 
when M2 — \fJ.\. 

II We will see in the next Section that the definition of the higgsino-like neutrahno does matter if 
M2 — |/i| or —Ml ~ |/i|. The definition we employ here maintains perturbative stability also in these 
scenarios with strong mixing. 
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In schemes la and lb this inversion can be done analytically [TO]. Here both chargino 
masses are inputs. The chargino mass matrix can easily be diagonalized analytically; this 
can be inverted to derive analytical expressions for M2 and /i in terms of the chargino 
masses (and tan/3, which we assume to be known independently). In general there are four 
solutions for M2 and fi. To begin with, the eigenvalue equations (of M^M^'^ or M'^^M^) 
are symmetric under M2 -H- fi. This degeneracy can be lifted once we know whether the 
lighter or the heavier chargino is (more) wino like; this determines whether M2 is smaller 
or larger than |/i|. Notice that even in this scheme knowledge of the qualitative properties 
of the chargino states is required to complete the inversion. The second degeneracy occurs 
since the sign of /i has not been fixed; recall that we define M2 to be real and positive. 
This degeneracy can be lifted by measuring some chargino coupling or cross section [23] 
in addition to the masses, since the two solutions will have different mixing matrices 
U and V. Alternatively one can use information from the neutralino sector to lift this 
degeneracy, since the two solutions will lead to somewhat different neutralino spectra. To 
this end, one can simply try both solutions, and see which one more accurately reproduces 
the measured neutralino spectrum. 

Having determined M2 and fi from the chargino sector, the remaining parameter Mi 
can be determined analytically from the neutralino sector [10]. Again, there is a two- 
fold degeneracy, having to do with the sign of Mi. This degeneracy would not occur if 
we were able to determine the eigenvalue of the neutralino mass matrix (jH]) rather than 
only its absolute value, which is the physical mass. In fact, the relative sign between 
two eigenvalues, with indices i and j, of this matrix is physical [5]: if the relative sign is 
positive, the Zx°X? coupling will be purely axial vector, whereas for negative relative sign 
it will be purely vector. This can be determined through e~^e~ — )■ x^Xp where a vector 
(axial vector) coupling leads to an S* — (P— )wave behavior of the cross section, i.e. to 
a suppression of the cross section by one (three) factors of the neutralino center-of-mass 
three-momentum. 

In the remaining schemes an analytical inversion from input masses to input param- 
eters is not practical. A numerical solution of the inversion equations is straightforward, 
though. Again these equations will have several solutions, which differ by the signs of 
Ml and/or /x, and/or by the ordering of the absolute values of Mi, M2 and /x. The latter 
degeneracy can again be lifted if information about the characters of the states is avail- 
able. For example, |Mi| will generally be smaller than M2 if the bino-like neutralino is 
lighter than the wino-like one, and so on. The sign issue can be resolved as in scheme 1. 

Note that in this manner one determines Mi, M2 and /i in an on-shell scheme, as 
required for our calculation. In order to relate these to DR parameters, one needs infor- 
mation about the entire superparticle and Higgs spectrum 



3.3 Combination with Spectrum Calculators 

Predictions for the masses of superparticles and Higgs bosons for given values of the 
relevant parameters of the MSSM Lagrangian are nowadays usually obtained with the 
help of publicly available programs like ISASUSY [25], SuSpect [22], SOFTSUSY [27] 
and SPHENO |28j. As per the SUSY Les Houches Accords (SLHA) [291 EQ] and the 
Supersymmetric Parameter Analysis (SPA) convention [31], values for parameters in the 
Lagrangian are given in the DR scheme, not in the on-shell scheme. This is true also for 
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the "weak-scale" output parameters Mi, M2 and /i of the spectrum calculators (which are 
typically defined at a SUSY breaking scale rather than the weak scale). These parameters 
should not be used as inputs in the mass matrices ([3]) and ([8]), since our input parameters 
are on-shell parameters; that is, in the limit of no mixing, where all differences between 
I Ml I, M2 and |/x| are much bigger than M^, these parameters become physical (on-shell) 
masses in our scheme, which is not the case for the running (DR) parameters. 

Fortunately modern spectrum calculators already perform the transformation from 
DR to on-shell parameters internally, in the limit of no mixing in the neutralino and 
chargino sector. Hence one can identify our on-shell input parameters Mi, M2 and —fi 
with whatever the spectrum calculators input as the (1,1), (2,2) and (3,4) entries of 
the neutralino mass matrix, respectively. Alternatively one can perform the conversion 
oneself, using the results of ref. [24]. 

However, spectrum calculators generally do not use on-shell W and Z masses in the 
off-diagonal elements of the chargino and neutralino mass matrices. In fact, physically 
there is no reason why on-shell masses should be preferable here. We have used them 
since the on-shell renormalization of the electroweak sector is very well established, and 
since we do not expect large radiative corrections from these sources anyway. Similarly, 
spectrum calculators do not enforce cos 9w = Mw/Mz, which in our scheme is necessary 
since otherwise some divergencies remain. The off-diagonal entries of the chargino and 
neutralino mass matrices [other than the (3,4) and (4,3) entries of the latter] should 
therefore not be copied from the spectrum calculators, but have to be entered "by hand" 
following Eqs.([3]) and ([8]). For this reason, the mixing matrices U, V and A^ should not 
be taken from the spectrum calculators, either; rather, they should be calculated by 
explicitly diagonalizing the chargino and mass matrices, constructed as per the above 
description, e.g. by using some standard matrix diagonalization routines. 

4 Numerical Results 

In this section we numerically compute the physical masses of the charginos and neu- 
tralinos, given by Eqs. flT5]) and fITB]) . We have used FeynArts [321 E3], FormCalc [21] 
and LoopTools [32] for the computations of the relevant two-point functions. Feynman 
gauge has been used throughout. For regularization we have used the constrained differ- 
ential renormalization method [36]. At one loop level this method has been proved to be 
equivalent [35] to regularization by the dimensional reduction method [37] . All numerical 
examples presented below employ the following parameters: 

A = 400 GeV, MgusY = 450 GeV, tan/3 = 5, M^ = 250 GeV. 

For simplicity we took the A parameters and the supersymmetry breaking soft sfermion 
masses Msusy to be universal as well as flavor conserving. 

In the flrst Subsection we survey the parameter space by scanning /i or Mi. We will 
see that some counterterms can become singular, leading to ill-deflned one-loop masses 
at the pole and, in some cases, to very large corrections near the pole. This issue is 
discussed in depth in the second Subsection. The flnal Subsection contains numerical 
comparisons between schemes for a couple of benchmark points. 
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Figure 1: Radiative corrections (in %) to all the neutralino masses against fi. The masses 
of two charginos and one neutralino have been used as inputs. Frame (a) is for scheme 
la, where the mass of the lightest neutralino is used as input, whereas frame (b) is for 
scheme lb, where instead the mass of the most bino-like neutralino is used as input. 



4.1 Surveys of Parameter Space 

In the following figures we show the differences between one-loop and tree-level masses 
in percent, defined by 



A? 



A; 



1— loop, OS tree 

Xi 
1— loop, OS tree 

^tree 



X 100, iG {1,2,3,4} 



X 100, ie{l,2} 



(22) 



Of course, these differences vanish for the three input states. We have fixed Mi = 
89.29 GeV, Ms = 369.78 GeV in all the plots against /x, and Ms = 200 GeV, /i = 300 GeV 
in all the plots against Mi. 

In Figure [1] we plot the corrections to the neutralino masses in scheme 1 against /i, 
with the left (right) frame corresponding to scheme la (lb). Both frames show singular- 
ities at /i = Ms. The reason is that, in these schemes, the expressions for 6M2 and Sn 
are undefined at Ms = |/i|, because 6M2 and 6fi are proportional to (M| — /^^) [T^[T3]. 
Hence the renormalized masses are undefined if Ms = |/i|. Moreover, some corrections 
remain relatively large as long as Ms and |/x| are close to each other. In particular, the 
mass of X3, which is mostly higgsino-like in this region, receives a large correction, mostly 
from 6fi. This pole is rather broad, e.g. A3 can exceed 20% if ||/i| — Ms| < 5 GeV. 

Figs, d] also show discontinuities at /x ~ 131.50 GeV. Here a level crossing occurs, i.e. 
the absolute value of the negative eigenvalue of the neutralino mass matrix M'^, which 
used to give m^o, decreases past a positive eigenvalue of M^ to give fn^o. However, if we 
track the corrections corresponding to the negative eigenvalue, there is no discontinuity. 
A similar situation is encountered for /i ~ —100 GeV, which is not shown in our plots. 
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Figure 2: Radiative corrections (in %) to all the neutralino masses against Mi. Frame 
(a) is for scheme la, and frame (b) is for scheme lb. Note that the two frames employ 
different scales along the ?/— axis. 



With this set of parameters, schemes la and lb perform equally well, i.e. the correc- 
tions have similar magnitudes. The results are different for |/i| < Mi, when the lightest 
neutralino is no longer bino-like. This leads to another discontinuity in scheme lb at 
yU ^ 110 GeV: for smaller (larger) fi the most bino-like neutralino is Xs (Xi)- When 
approaching this point from above, Ag therefore switches from a positive value to zero, 
while A5 becomes non-vanishing (and negative, for this choice of parameters) . Moreover, 
(5 Ml becomes large in scheme la if \n\ <^ Mi, as anticipated in the discussion in Sub- 
sec J2.3j see also Table [TJ This leads to a sizable correction to the mass of the bino-like 
neutralino. 

This does not cause serious instability in Fig. [H since the small value of Mi chosen 
implies that \Mi — fi\ cannot exceed Mz, i.e. Xi retains a significant bino component 
for all values of /i shown. Note also that LEP searches for xfxi and X1X2 production 
anyway exclude the region fi < 120 GeV. The problematic behavior of scheme la for 
|yu| < Ml — Mz becomes more evident in Fig. |2l where we vary Mi instead of /x. As Mi 
grows large, the bino component of the lightest neutralino decreases, Nu oc Mz/{Mi—fj,). 
Hence 6M1 becomes large in scheme la, giving a large correction to the mass of the bino- 
like neutralino, which is xt fo^^ ^1 > ^2 in this figure. Note that in anomaly-mediated 
supersymmetry breaking one expects Mi ~ 2.8M2 [3]. As a consequence \Mi\ — \fi\ > Mz 
can easily be realized in such a scenario. 

Fig. [2] shows some discontinuities at Mi ~ 335 GeV. Just below this value X2 is the 
most bino-like neutralino, but at Mi = 335 GeV Xa becomes more bino-like. Shortly 
thereafter, at Mi ~ 340 GeV, there is a level crossing between X2 ^^d x^; the corre- 
sponding eigenvalues of the neutralino mass matrix have opposite sign. Note that all 
corrections remain below 1% in magnitude near these discontinuities. 

In Fig. 121 we demonstrate that scheme 2b also performs equally well with our choice of 
parameters. However, the counterterms again show a singularity at M2 — l/^l (although 
not exactly at M2 = |/i|), where the denominator D of the counterterms, given in Eq. (l26|) . 
vanishes. This again leads to large corrections to the mass of the higgsino-like neutralino 
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Figure 3: Radiative corrections (in %) to (a) the neutralino masses and (b) the chargino 
masses against /i. Results are for scheme 2b, where the masses of the higgsino-hke 
chargino and of the wino- and bino-hke neutrahnos are used as inputs. 
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Figure 4: Radiative corrections (in %) to (a) the neutrahno masses and (b) the chargino 
masses against /x. Results are for scheme 2c, where the masses of the wino-like chargino 
and of the bino- and higgsino-like neutralinos are used as inputs. 

Xg in the vicinity of the singularity, very similar to both versions of scheme 1. 

In contrast, we find that scheme 2a (not shown) performs poorly whenever M2 and 
|/i| are very different, as anticipated in the discussion of Table [1] In fact, numerically 
it is considerably worse than scheme la. In the latter only SMi is occasionally poorly 
determined; the corresponding corrections are determined by the f/(l)y coupling. In 
scheme 2a, S^ or SM2 can be poorly determined, which receive corrections proportional 
to the larger SU{2) coupling. As a result the region of parameter space with large 
corrections is not only larger in scheme 2a than in scheme la, the size of the correction 
in these regions also tends to be larger. 

Scheme 2c does not have any singularity near M2 = \fi\. In Fig. [Hit produces stable 
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Figure 5: Radiative corrections (in %) to (a) the neutralino masses and (b) the chargino 
masses against Mi. Results are for scheme 2c, where the masses of the wino-hke chargino, 
bino- and higgsino-hke neutrahnos are used as inputs. 
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Figure 6: Radiative corrections (in %) to (a) the neutrahno masses and (b) the chargino 
masses against fi. Results are for scheme 3a, where the masses of the lightest three 
neutrahnos are used as inputs. 

and small corrections to all three chargino and neutralino masses that are not inputs for 
the entire range of /i. The curves in this figure show several discontinuities, but these are 
harmless, since they do not endanger perturbative stability. In particular, in frame (a) 
we observe the same discontinuity as in Fig. [H since the states Xs and X2 cross each other 
at this point. The discontinuities at yU ~ 110 GeV have the same origin as in scheme lb, 
see the discussion of Fig. [U^b). Finally, at /i ~ 370 GeV (i.e. M2 — /x), the wino-like 
state changes from xt ^o xt as /i increases. This affects 6M2 significantly, giving rise to 
(harmless) discontinuities in both frames of Fig. |H 

In Fig. Owe illustrate results for scheme 2c as function of Mi. The discontinuities 
in the corrections to neutralino masses near Mi = 335 GeV have the same origin as in 
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Figure 7: Radiative corrections (in %) to (a) the neutralino masses and (b) the chargino 
masses against /i. Results are for scheme 3b, where the masses of the bino-, wino- and 
higgsino-hke neutrahnos are used as inputs. 
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Figure 8: Radiative corrections (in %) to (a) the neutrahno masses and (b) the chargino 
masses against Mi. Results are for scheme 3b, where the masses of the bino-, wino- and 
higgsino-like neutrahnos are used as inputs. 

scheme lb, see Fig. Et^b). All corrections are well behaved for |Mi| ~ |yu|. 

Fig. [6] demonstrates that scheme 3a also behaves poorly. For /i ^ 200 GeV all three 
input states, which are the lightest three neutrahnos in this scheme, have small wino 
component. As a consequence the 6M2 is poorly determined, leading to huge corrections 
to the masses of the wino-like states Xa and xt^ largely from SM2. The scheme also shows 
instability around Mi = /i, giving very large corrections to the masses of both charginos 
and the wino-like neutralino. 

Scheme 3b does not have any singularity near M2 = \n\. In Fig. [71 like scheme 2c, it 
also produces stable and small corrections to all three chargino and neutralino masses that 
are not inputs for the entire range of /i. In frame (a) we observe the same discontinuity 
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as in Fig. [H since the states Xs ^i-nd X2 cross each other at this point. The discontinuities 
at /i ~ 110 GeV have the same origin as in scheme lb, see the discussion of Fig. [It^b). 
The harmless discontinuities at /x ~ 370 GeV (i.e. M2 — fi), are essentially the same as 
in Fig. m 

In Fig. [H] we show results for scheme 3b as function of Mi. The discontinuities in 
the corrections to neutralino masses near Mi = 335 GeV have the same origin as in 
scheme lb, see Fig. Mjo)- In the case at hand this also leads to a small discontinuity 
in A]''. All corrections are well behaved for \Mi\ ~ |/i|, but there is some instability 
around Mi = 178.14 GeV. Here the denominator D, given in Eq.( l33|l . of the counterterms 
vanishes. As a result, the neutralino and chargino masses become ill-defined at this point. 
The effect of this pole is most visible in the mass of the lightest chargino which, being 
wino-like, receives dominant correction from 6M2. However, even this corrections drops 
to ~ 2% when Mi is just 1 GeV away from the pole, i.e. this pole is much narrower than 
the poles at M2 = |/i| in schemes la and lb, and the poles at M2 — |/i| in schemes 2a and 
2b. These poles are discussed in more detail in the following Subsection. 

4.2 Singularities in the Counterterms 
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Figure 9: Contours of constant radiative corrections to (a) the higgsino-like neutralino 
mass around M2 ~ /^ in scheme lb, (b) the wino-like chargino mass around Mi ~ M2 in 
scheme 3b have been plotted against /x and Mi respectively. Only corrections exceeding 
10% have been shown. The patch in (b), at Mi ~ 300 GeV, is actually at Mi ~ —300 
GeV, i.e. it is plotted against —Mi. 



In Fig. in] the occurrence of poles, i.e. true divergences in the "finite" parts of the 
counterterms (as opposed to harmless discontinuities due to level crossings and the like), 
is analyzed further. These are not due to divergent loop integrals, but due to zeroes of 
the denominators in the expressions for the counterterms [Eqs. fj26|) and fl33|) . and the 
corresponding expression for scheme 1]. We have assumed Mi = 180 GeV and fi = 300 
GeV in Figs. IHti^ andlHt", respectively. Here we show the maximal (relative) corrections. In 
scheme lb these occur for the higgsino-like neutralino, i.e. Fig. [9^ shows Ag (for fi ^ 168 
GeV as well as 168 GeV < /i < 174 GeV and M2 > /x) or A|] (for /i > 174 GeV as well as 
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168 GeV ^ yU ^ 174 GeV and M2 < jj)\j In scheme 3b the mass of the wino-hke chargino 
instead receives the largest correction, i.e. Fig. [9b, shows A^ (for M2 > /i) or A^ (for 
M2 < n). Finally, scheme 2b shows instability both at M2 — |/u| and at Mi ~ M2. 

As already mentioned, it is easy to see that 6M2 and 6fi diverge in both versions of 
scheme 1 at M2 = |/i| [HlIIS]. In schemes 2 and 3 the expressions for the denominators of 
the counterterms, given in Eqs. (l26|) and (133|) . respectively, are too complicated to allow 
the derivation of analytical expressions for the pole positions in terms of the input pa- 
rameters Ml, M2, n and tan /3. However, it is easy to see that in the limit of small mixing 
only one term in each of these expressions is large, leading to non-vanishing (and sizable) 
denominators. Cancellations leading to poles can only occur in the presence of large mix- 
ing, which requires at least one of the differences ||Mi| — M2I , ||Mi| — \fi\\ , IM2 — |Ai|| 
to be relatively small. Note that this is a necessary condition, not a sufficient one. 

In fact, Eqs. fl26|l and f l33|) show that singularities can only occur if at least two of 
the input states are strongly mixed. This explains the absence of poles at M2 ^ |/i| 
in schemes 2c and 3b. Here the most higgsino-like neutralino is chosen as input state. 
Note that even at M2 = |yu| only one of the two higgsino-like neutralinos mixes strongly 
with the wino; the other state corresponds to an eigenvalue with opposite sign (~ — |/x|), 
and shows little mixing. Small mixing means large higgsino component, i.e. this weakly 
mixed state is picked up by our algorithm as the most higgsino-like one. At M2 = \n\ (and 
sufficiently large ||Mi| — M2I) the input states in these schemes are thus a relatively pure 
bino (neutralino), a relatively pure higgsino (neutralino), and a strongly mixed wino- 
higgsino state (a chargino in scheme 2c, a neutralino in scheme 3b). The fact that two 
input states show little mixing is sufficient to single out one term in eqs. (l26p and fl33l) as 
dominant one, making significant cancellations impossible. A similar argument explains 
why none of our schemes has instabilities for |Mi| — |/x| (unless one also has M2 — |Ai|; 
see below). 

Scheme 3b shows a singularity at Mi ~ M2 where two input neutralino states are 
strong mixtures of the bino and wino current states, leading to a cancellation between 
the last two terms in Eq. (l33l) with all other terms being small. Note that the pole occurs 
only when |Mi - |/i|| > IM2 - \fi\\^ We did not find any pole if |Mi - /i| < 80 GeV. 
In this interval one of the bino-wino mixed states (usually the more wino-like one) also 
has a substantial higgsino component, and correspondingly reduced gaugino components. 
Thus there are no longer two input states with very similar bino and wino components. 

This pole at Mi ~ M2 is very narrow. This is partly due to cancellations of the 
numerators of the expressions for the counterterms, which can be understood as follows. 
Near the pole, indices i and j denote strong mixtures between bino and wino states, 
and k denotes a higgsino state with little mixing (unless ||yu| — M2I is also small). In 
that case only two terms in each of Eqs. (l34l) - (l36l) are not suppressed. These terms are 
antisymmetric in i -H- j, hence strong mixing also induces strong cancellations; however, 
these cancellations are not perfect, since mixing pushes the masses m^o and m^o apart. 



*Iii the region 168 GeV < /i < 174 GeV the higgsino-hke neutrahno trades places with a photino-hke 
state. Since these states have eigenvalues of the neutralino mass matrix with different sign, this cross- 
over of states is not associated with strong additional mixing, and thus does not cause any feature in 
Fig. Hi. 

tPor exact equality, Mi — M2, one eigenstate is an exact photino with mass M2, whose bino- 
component is given by cos 0w ~ 0.88, i.e. with relatively weak bino-wino mixing. 
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The finite parts of tlie counterterms tliemselves can tlierefore become quite large, often 
exceeding the tree-level values of these parameters. Even then corrections to physical 
masses need not be large, due to cancellations in Eqs. fllSI) and flT6|) . By construction 
these cancellations are perfect for the three input states. If an output state is of similar 
nature as one of the input states, there will then also be strong cancellations for the mass 
of this output state. For example, in scheme lb with M2 — |/i|, the strong wino-higgsino 
mixtures in the neutralino sector have similar wino and higgsino components as the two 
input chargino states, leading to strong cancellations. The largest correction then results 
for the (almost) unmixed higgsino state in the neutralino sector, as shown in Fig.[9K. The 
same argument explains why in the "generic" situation with relatively small mixing, in 
all our b and c schemes the mass of the most strongly mixed state will receive the largest 
correction, since it is most different from all the input states, which are the least mixed 
ones in these schemes. 

In scheme 3b with Mi ~ M2, both strongly mixed gaugino-like neutralinos, which 
are the cause of the pole, are input states, and hence have perfect cancellations. These 
states mix relatively weakly with the higgsino-like states. The counterterm 5yU, which 
diverges at the pole, comes with coefficient —2Nl^N^^ in Eq. (IT6|) . In the limit of vanishing 
higgsino-gaugino mixing, this is the same for the two higgsino-like states, one having 
iVfc3 ~ Nk4 — i/V^ (purely imaginary), while the other had N^^ ~ —N^^ = l/\/2 (purely 
real, but with opposite signs). This explains why the correction to the mass of the single 
output state in the neutralino sector remains very small near this pole, see Fig. El 

One of the two chargino states is also higgsino-like, and the cancellation will go 
through. However, the second chargino is wino-like. Since there is no wino-like neu- 
tralino, only strongly mixed ones, the cancellation in the equation for the mass of this 
wino-like chargino is quite incomplete. This explains why this state receives the largest 
correction, as shown in Fig. [Hb- 

The cancellations we discussed so far work (to some extent) near all poles, and help 
to confine their effects to quite narrow strips of parameter space, as shown in Fig. [91 
The fact that the instability region of scheme 3b at Mi ~ M2 is much narrower than 
that of scheme lb (or 2b) at M2 — |/i| is due to the fact that wino-bino mixing, which 
causes the former, only occurs at second order in an expansion of the neutralino masses 
and eigenstates in powers of Mz, while wino-higgsino mixing, which causes the latter, 
already occurs at first order; wino-bino mixing is also suppressed by an extra factor of 
sin 6w- As a result, significant wino-bino mixing occurs over a much smaller region of 
parameter space than large wino-higgsino mixing. 

Finally, if —Mi ~ M2 — \fi\ there are two bino-higgsino mixed states and two wino- 
higgsino mixed states in the neutralino sector. Thus scheme 3b also suffers from instability 
around this point because two input states are strongly mixed|i| This explains the patch 
in Fig. [9]d at Mi ~ —300 GeV. Note that here the "finite" parts of the counterterms 
diverge only in a very small region in the (Mi, M2) plane, extending for about 1 MeV in 
M2 and about 300 MeV in Mi; this instability is thus qualitatively different from the one 
at Ml ~ M2, where a pole occurs along a line that is unbounded towards large M2la 



■l-For tan f3 — I there is a pure higgsino state in the neutrahno seetor. The eorresponding eigenvalue is 
given by — /i. Thus, even with —Mi ~ M2 — |/i|, as tan /3 approaches 1, there is only little bino-higgsino 
(wino-higgsino) mixing if /i > (/x < 0). 

§In this patch the two neutralino states that remain after identifying the most bino- and wino-like 
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The singularity at Mi ~ M2 does not exist in scheme 2c, since in this scheme again 
only one of the three input states, the most bino-like neutralino, is strongly mixed in this 
region of parameter space. The input chargino state is mostly wino-like, unless one in 
addition has M2 ~ |//|, and the second input neutralino is again mostly higgsino (even if 
M2 — |/i|, as explained above). 

For fixed /i and tan/3 Scheme 2c can therefore at worst become singular in a very 
narrow patch in the (Mi, M2) plane, in the region —Mi ~ M2 — |/i|. We saw above that 
this leads to two bino-higgsino mixed states and two wino-higgsino mixed states in the 
neutralino sector. The chargino sector also contains two strongly mixed (wino-higgsino) 
states. Note that the denominator D of Eq. 026p is not a continuous function of the 
input parameters, since the indices i,j, k defining the input states take different values in 
different regions of parameter space. This makes it rather difficult to systematically search 
for zeroes of this denominator. A straightforward numerical scan of Mi and M2 with step 
size of 0.01 GeV did not find any evidence for a zero of the denominator. We did find 
points with large (finite parts of) SMi, 6M2 and 5/i in a very narrow region of parameter 
space, similar to the patch in Fig. [Hb- However, even in that region the renormalized 
masses remain well behaved in scheme 2c, the maximal correction not exceeding 4%. 
This is due to the two types of cancellation discussed above. Note that all our other 
schemes lead to considerable larger corrections around this point. We thus conclude that 
scheme 2c is the most stable of the seven schemes we investigated. 

4.3 Comparison Between Schemes 

It is generally expected that the physical masses calculated using two different renor- 
malization schemes at one-loop level can differ only by terms of two-loop order. In the 
case at hand this is true whenever the schemes are well behaved. We demonstrate this 
explicitly by comparing results for schemes lb and 3b. We then show another example 
where this no longer holds, since one of the schemes produces large corrections close to 
one of the singularities discussed in the previous Subsection. We do not discuss schemes 
la, 2a and 3a in this Subsection, since we already saw that they perform poorly in large 
regions of parameter space, see Table [TJ 

We use the following "tree level" parameters to calculate the on-shell renormalized 
masses in scheme (lb): 

Set A: Ml = 90 GeV, M2 = 200 GeV, /i = 300 GeV, tan /3 = 5. 

The masses thus obtained are shown in column 3 of Table [2l We then take m^o , m^o and 
m-j^o, thus obtained, as input masses for scheme 3b. Note that these are respectively the 
masses of the most bino-, wino- and higgsino-like neutralino. From these three masses, 
we compute numerically Mi, M2 and /i: 

Set B: Ml = 89.9697 GeV, M2 = 199.7935 GeV, /i = 301.0052 GeV, tan/3 = 5. 



neutralinos have very similar higgsino components. The instability can be avoided by simply using the 
other higgsino-like state as input state, even though by our definition it is slightly less higgsino-like. 
Note also that disentangling the states experimentally will be very difficult in this region, again due to 
the strong mixing. 
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We then calculate the one-loop masses in scheme 3b. The resulting spectrum is shown 
in column 4 of Table [2j We also list in parentheses the one-loop masses calculated using 
parameter set A within scheme 3b. As shown in column 5, the differences between the 
two one-loop predictions, using scheme lb with Set A and scheme 3b with set B, are 
much smaller than the differences between tree-level and one-loop corrections (except 
for the input masses, of course, where the latter differences vanish). This is a non-trivial 
check of our calculation. 

Note that using the same numerical values for Mi, M2 and /z in both schemes leads 
to differences between one-loop masses that are themselves of one-loop order, as shown 
by the numbers in parentheses in column 5 of Table [2J Only after fixing three physical 
masses to be the same in both schemes do the predictions for all other masses become 
nearly identical. The reason is that Mi, M2 and /i are themselves scheme dependent 
quantities. 



Particle 


mf" [GeV] 


1- 


-loop, OS l-QgYj 


dniy^ [GeV] 




[set A] 


scheme lb 

[set A] 


scheme 3b 
[set B (set A)] 


col. (3) -col. (4) 


xt 


170.4007 


170.4007 


170.4030 (170.3714) 


-0.0023 (0.0293) 


xt 


337.5026 


337.5026 


337.5169 (336.7472) 


-0.0143 (0.7554) 


x\ 


84.8425 


84.8425 


84.8425 (84.8425) 


0.0 (0.0) 


xl 


172.1478 


172.1545 


172.1545 (172.1478) 


0.0 (0.0067) 


x^ 


305.4830 


306.4778 


306.4778 (305.4830) 


0.0 (0.9948) 


xl 


338.4927 


339.6138 


339.6313 (338.8642) 


-0.0175 (0.7496) 



Table 2: Masses of the charginos and neutralinos. The second column lists tree-level 
predictions for input set A. Columns 3 and 4 list one-loop predictions for schemes lb 
and 3b; in the latter case, the first entry is for input set B, which has been chosen such 
that Tn^Q, m^o and m^o have the same values as in scheme lb (third column), whereas 
the numbers in parentheses use input set A. The last column is the difference between 
the third and fourth columns. 

However, as mentioned before, scheme lb is undefined when M2 = \fi\. We also saw 
that for M2 — |/i| the finite part of the counterterm 6fi remains large, so that the mass of 
a higgsino-like neutralino receives a large correction. This pole is absent in scheme 3b. 
As a consequence, for the same set of physical masses, the one-loop masses in these two 
schemes can be quite different. In Table [3] we demonstrate this by explicit computation. 

We use the following tree level parameters to calculate the on-shell renormalized 
masses in scheme lb: 

Set C: Ml = 90 GeV, M2 = 350 GeV, /i = 355 GeV, tan/3 = 5. 



The resulting physical masses are given in column 3 of table 
correction (> 50 GeV) to m^o|j As before, we then take 



m^o, 171^2 and m^o 

Xi' Xo X3 



Note the very large 
thus obtained, 



^In case of rriyO and m^o large finite corrections from 6M2 and (5/i largely cancel, similar to the exact 
cancellation required to keep both chargino masses unchanged. X3 is ^ nearly pure higgsino; its mass is 
thus not sensitive to 5M2, so that no cancellation can take place. 
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as input masses for scheme 3b. These are respectively the masses of the bino-, wino- 
and higgsino-hke neutrahno. From these three masses, we compute numerically Mi, M2 
and fi: 

Set D: Ml = 89.2056 GeV, M2 = 319.1479 GeV, /x = 418.1529 GeV, tan/3 = 5. 

Note that fi had to be increased by more than 60 GeV in order to reproduce the one-loop 
prediction for m^o from scheme lb. This increase of /i would also have led to sizable 
increases of m-+ and m^o, had we not simultaneously reduced M2 by more than 30 GeV. 
Only Ml retains approximately the same value in both Sets. 



Particle 


m*/^*' [GeVl 


m^- 


-loop, OS FQgyl 


(5m^ [GeV] 




[set C] 


scheme lb 

[set C] 


scheme 3b 

[set D (set C)] 


col. (3) -col. (4) 


xt 


288.3284 


288.3284 


288.0337 (287.9178) 


0.2947 (0.4106) 


xt 


422.2986 


422.2986 


452.4840 (417.3805) 


-30.1854 (4.9181) 


X? 


86.3675 


86.3675 


86.3675 (86.3675) 


0.0 (0.0) 


xl 


289.8939 


288.7260 


288.7260 (289.8939) 


0.0 (-1.1679) 


xl 


359.1238 


421.9883 


421.9883 (359.1238) 


0.0 (62.8645) 


xl 


422.8623 


424.2374 


454.7079 (418.8926) 


-30.4705 (5.3448) 



Table 3: Masses of neutralinos and charginos. The meaning of the various entries is as 
in Table El except that parameter sets A and B have been replaced by sets C and D, 
respectively. 

The resulting one-loop masses in scheme 3b are listed in column 4 of Table [31 In 
parentheses we also give the one-loop masses calculated using parameter set C using 
scheme 3b. Column 5 shows that the differences between the one-loop predictions for 
m^+ and m^o are very large. This is due to the large differences in the input parameters 
between Set C and Set D, which in turn was enforced by the large correction to m^o in 
scheme lb. 

We emphasize that in the case at hand, the predictions of scheme lb are clearly not 
reliable, due to the vicinity of the pole at M2 = /i which leads to a breakdown of the 
perturbative expansion in this scheme. In scheme 3b the differences between tree-level 
and one-loop masses do not exceed 1% for both set C and set D. Moreover, differences 
between physical masses are again very small if scheme 3b is compared to scheme 2c, 
which also has no pole at M2 ~ /i. 



5 Conclusion 

In this paper we have analyzed variants of the on-shell renormalization of the chargino 
and neutralino sector of the MSSM. There are six physical particles with in general 
different masses in this sector, but only three free parameters whose counterterms are 
determined from this sector. Hence only three masses can be chosen as "input masses", 
which by definition are unchanged to all orders in perturbation theory; the remaining 
three masses will receive quantum corrections, which we computed to one-loop order. 
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We introduced seven different versions of tlie on-sfiell scfieme in Sec. 2. Tliese scliemes 
differ by tlie states wfiose masses are used as inputs. General considerations led us to 
expect that all schemes where the indices of the three input states are fixed a priori can 
show instability over large regions of parameter space, as listed in Table [T] for schemes 
la, 2a and 3a; we confirmed this by explicit calculations in Subsec. 4.1. In these regions 
at least one of three counterterms 5 Mi, SM2 and 6^ is poorly determined, since none of 
the input masses is very sensitive to it. These large regions of instability can be avoided 
by choosing the masses one bino-like, one wino-like and one higgsino-like state as input; 
the definition of these states, and other practical matters, have been discussed in Sec. 3. 
This ensures that there is at least one input mass that is sensitive to each counterterm. 

Unfortunately this is not sufficient to avoid all regions of instability. In particular, 
we see that all schemes where the mass of the higgsino-like chargino is used as an input 
become unstable for M2 — \n\, because finite parts of the counterterms 6M2 and 6fi 
become very large, leading to a very large correction to the mass of the (most) higgsino- 
like neutralino. This instability is caused by numerical cancellations in the denominators 
of the expressions for the counterterms, which can only occur if at least two input states 
are strongly mixed states. We saw in Subsecs. 4.2 and 4.3 that this instability can lead 
to a collapse of the perturbative expansion over a significant region of parameter space. 
It can be avoided if instead the mass of the higgsino-like neutralino is used as input; this 
state does not mix strongly even if M2 — \n\ or \Mi\ — |/x|. 

Finally, we found another similar pole, with however in general much smaller numerical 
values of the corrections somewhat away from the pole, if the masses of both the bino- 
and wino-like neutralino are used as inputs. This instability is due to bino-wino mixing, 
which occurs at Mi ~ M2. It can be avoided by chosing the mass of the wino-like 
chargino, rather than the wino-like neutralino, as input, since the charginos are obviously 
not affected by bino-wino mixing. 

We conclude that the perturbatively most stable version of the on-shell scheme is 
the one where the masses of the wino-like chargino, and of the bino- and higgsino-like 
neutralinos are used as inputs; the scheme where the masses of the bino-, wino- and 
higgsino-like neutralinos are used as input performs as well over most, but not all, of 
parameter space. 

We emphasize that a poor choice of renormalization scheme can lead to unphysically 
large corrections not only to some masses, but also to cross sections. For example, the 
exchange of heavier neutralino states can play an important role in the cross section for 
the annihilation of a pair of the lightest neutralinos into final states containing gauge 
and/or Higgs bosons [M]; the total neutralino annihilation cross section determines their 
present relic density if they are stable. The proper choice of renormalization scheme 
should therefore be important for a variety of one-loop calculations in the MSSM that 
have been, and are being, performed. 
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Appendix A 

Here we give expressions for 5 Mi , 5M2 and 5^ in terms of the 5my. of Eq. (TT2|) as well as 
the counterterms to the off-diagonal entries of the chargino and neutralino mass matri- 
ces, which can be derived from Eqs. (IT7|) and (TT8l) . These expressions have been obtained 
by solving equations (120|) . i.e. they are valid for all versions of scheme 2, with differ- 
ent variants corresponding to different choices of the indices i,j, A; (see Tabled]). The 
counterterms can be written as 

5Mi = ^, (23) 

SM, = ^, (24) 

S, - §, (25) 

where 

D = 2n:,n:,n;Iu;,v:, - 2n:^ n;,n;,u;,v:, + n:^n;Iu;,v:, - NtNpi^v,,, (26) 

N, = - {{5m^. - 25min:,n:, - 26m^,n:,n:, - 25m^,n:,n:, 



-25Mli,N:^N:,)N*l + N:l{-5m^. + 25M^,N';,N*^ + 25M:i^N;,N;, 



(27) 



Nm, = 26m^^N:,N*,Nf,-2Sm^.N:^N;,N;,-26M^,N:,N*,Nf,U;,V,\ 

+ 25m^,n:;n;,n;,u;,v:, - 25mi^n:,n:,n*Iui,v:^ - 6m^oN:;u;,v:, 

- 2sm^,n:,n:,n;Iu:,v:, - 25M^,N:,N:,Nf,u:,v:, - 25m^,n:,n:,n;Iu;,v:, 

+ 26m^,n*^n;,n;,u;,v:, + 26m^,n:; n;,n;,u:,v:, + 2Sm^,n*^n;,n;,ui,v:, 

+ 26m^,n:;n;,n;,u;,v:„ (28) 
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- 25M^,N:,N:,Nf,u:,v:, - 26m^,n:,n:,n;Iu;,v:, + 25m^,n:^n;,n;,u;,v:, 
+ 26m^,n*^n;,n*,u;,v:, + 26m^,n:^n;,n;,u;,v:, + 26m^,n*^n;,n;,u;,v,% 

- 26M^,,N*,N:,Nf,ui,v:, + 26m^,,n*^n*,n;MiV:2 + ^6 m^,n*,n*,n;,n;,ui,v,\ 
+ a6m^,n*,n:,n;,n*,u;,v:, - 26M^,N*,N*,N;iu;,v:, - 26M^,N*,N*,Nf,u;,v:, 

- a5m^,n:,n:,n;,n;,u;,v:, - 25m^^,n:^n;,n*Iui^v:, + 25M:i,N*iN*,N*,Ul2V:, 

+ 6m^.Nf,Ul,V:,. (29) 

Appendix B 

Here we give SMi, SM2 and Sfi for scheme 3, obtained from solving equations (!2T!) . 
Different variants of tliis scfieme correspond to different clioices of i,j,k (see Tabled]). 
The counterterms are given by 

SM^ = ^, (30) 

SM, = ^, (31) 

S, - §, (32) 



where 



D = {n*,n:,n;,n;,-n:,n;,n;,n*,-n:,n:,n;,n*,, + n:,n;,n;,n:, 

+NiN*lNl,Nl, - N:In;InI,NI,), (33) 

+ 26M^,N*,N:,N;^Nf, + 26M^,N:,N:,Nf,N*l + 26M^,N*,N*,Nf,N;l 

+ 26M^,N:,N*,N*,NII - 25M^,NiN*,N*,Nll - 26M^,N:^ N;,N;,NII 

- 26M^,N:;n;,N;,N;1 - 2SM^,N:^N;,N;,Nf, - 5m^.N:lNt, + 5m^.N;lNt, 

- 26M^,N:,N:,Nf,N:l - 25M^,N:,N:,Nf,Nf, - 25M^,N:,N:,Nf,Nf, 

- 26M^,N:,N:,N*In;1 + 26M^,N*;N;,N;,Nf, + 26M^,N*; N;,N;,N;1 

+ 26m^,n:;n;,n;,n;1 + 26m^,n:; n;,n;,n;1 + 26m^,n:^n;In;,n;, 

- 25M^,N:In*InI,NI, + 25Mli,NiN*lNl^Nl, - 26M^,N:^ Nf,N*,,N*,, 

+ 26m^,n:;n;In;,n;, - 26M^,N:^Nf,N;,N;, + 26m^,n:;n;In;,n; 



^2 , ,^2 , ,^ , ,^ „ , ,^2 , ,^2 , , ^ , 

- 26M^,N*;N;iN;,N;„ (34) 
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,2 



- 25m^,n:,n:,n;,n;,n:^ - 25m^,n:,n:,n;,n;,n:^ - 25m^,n:,n:,n;,n;,n:1 

- 26m^,n*,n:,n;,n;,n;1 - 6m^oN*,N*,N;l + 26m^,n*,n*,n;,n;,n;1 
+ 25m^,n*,n*,n;,n;,n;1 + 25m^,n*,n*,n;,n;,n;1 + 25m^,n:,n:,n;,n;,n;1 
+ 26m^,n*,n*,n*In;,n;, - 2Sm^,n:; n;,n;,n;,n;, + 2Sm^,n*,n:,n;In:,n:, 

- 26m-,n:;n;,n;,n;,n;, + 26m-,n:,n*,n;In;,n;, - 26m-,n:;n;,n;,n;,n;, 

- 25m^^n:,n:^n*Ini,ni, - 25m^,n:,n:^n*Ini,ni, - 25m^,n:,n:,n*Ini,ni, 

- 26m^,n*,n:,n;In;,n;, + 26m^,n*^n;,n;,n;,n;, + 26m^,n:;n;,n;,n;,n;, 
+ 26m^,n*;n;,n;,n;,n;, + 26m^,n*;n*,n*,n;,n;„ (35) 

and, 

Nm, = Sm^oN:,N:,Nf, - dm^oN*^ N;,N;, - dm^^Nl^Nl^Nf, + 25M^,N:,N:,N;,N*,Nf^ 

+ 25Ml^,N:^N:,N*^N;,Ni, + 25M^,N:,N:,N;,N;,Ni, + 25Ml^,N:,N:,N;^N;,Nf^ 

+ 5m^.N;,N;,Ni, - 26M^,N:,N*,N*,N*,Nf, - 26M^,N:,N:,N*,N*,Nf, 

- 25M]^,n*,n:,n;,n;,n;1 - 2Sm^,n*,n*,n;,n;,n;1 - 25m^,n:,n:,n*Ini,ni^ 
+ 25M^^NiN*,N;,Ni,Ni, - 2SM^,N:,N:,Nf,N;,N;, + 2Sm^,n*;n;,n;,n;,n;, 

- 26M^,N*,N*,Nf,N;,N;, + 2Sm^,n:; n;,n;,n;,n;, - 2SM^,N*,N:,Nf,N;,N;, 

+ 26M^,NgN;,N;,Nl,Nl, + dm^.N^Nl^Nl, - 6m^oNf,N;,N;, 

+ 25M-,N*,N;,N*InI,NI, + 25M-,N:,N*,Nf,Nl,Nl, + 25M^,N;,N:,N*InI,NI, 

+ 26M^,N*^N:,N*InI^NI, - 26M^^N*lN*,N*,Nl,Nl, - 26M^,N*; N;,N;,N*,,N*,, 

- 26m^,n*;n;,n*,n;,n;,-26m^,n*^n*,n*,n;,n;,. (36) 



Appendix C 

We show by explicit calculation that using the xtx^Z couplings, it is possible to deter- 
mine which of the two charginos is more wino-like, or, equivalently, whether M2 is bigger 
or smaller than |/i|. The relevant terms we compare are 

a = i-Vl - Vj + sl^r + i-Ul - \ul + sl)\ V ^ G {1,2}, (37) 

which determine the sum of the squares of the left- and right-handed couplings. These 
terms contribute, for example, to the chargino pair production cross section through 
Z— boson exchange in the s— channel in e"^e^ or pp collisions. Recall that we use the 
convention M2 > and assume vanishing CP phases in the chargino mass matrix. How- 
ever, /i can be either positive or negative. The unitary mixing matrices U and V are real. 
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They satisfy 

f/2. = f/2. f/2 = f/2 

^' = ^' ^? = ^' (38) 

where {i,j} G {1, 2}. 

Recall that, for given tan/3, we can solve the equations for the physical masses m-+ 
for the parameters M2 and /i, with a four-fold ambiguity. Part of this ambiguity is due 
to the fact that the expressions for the physical masses are invariant under the exchange 
of M2 and /i. Here we show that the quantities Ci of Eq.f l37|l can be used to resolve this 
ambiguity. We show this separately for positive and negative yU. 

• case 1 : /i > 

Let a, b be the possible solutions for M2, fi] a,b > 0. For {M2, /x} = {a, b}, let U and 
V be the unitary mixing matrices satisfying Eq.(jl]). Then, for {M2, fi} = {b, a}, the 
corresponding matrices are given by VO and UO respectively, where 

o-{r,). ,30) 

Using relations fl38|) . it is straightforward to observe that in M^, a ■^b corresponds 
to Ci ^ C2. 

• case 2 : /i < 

Let a,b be the possible solutions for M2, \fi\; a,b > 0. For {M2,/i} = {a, —b}, let 

U and V be the unitary mixing matrices satisfying Eq.(jl]). Then, for {M2, n} = 

_ J' 

{b, —a}, the corresponding matrices are given by VO and UO respectively, where 

0=(_°ij). (40) 

Again, {a, —b} ^ {b, —a} corresponds to Ci ^ C2. 

We thus see that by comparing Ci with C2 it is possible in both cases to determine 
the ordering of M2 and \fi\. However, if M2 — \n\, it can be difficult to infer whether 
M2 > |/i| or vice versa. This is of some importance, since precisely in this region of 
parameter space schemes where the mass of the more higgsino-like chargino is used as 
input show an instability, as we saw in Sec. 4. 
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